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2n ■ Abstract 
G^ 
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. We present a general method for incorporating an external electromagnetic field 

, into descriptions of few-body systems whose strong interactions are described 

by integral equations. In particular, we address the case where the integral 
equations are those of quantum field theory and effectively sum up an infinite 
number of Feynman diagrams. The method involves the idea of gauging the 
integral equations themselves, and results in electromagnetic amplitudes where 
, an external photon is effectively coupled to every part of every strong interaction 

I diagram in the model. Current conservation is therefore implemented in the 

Q I way prescribed by quantum field theory. We apply our gauging procedure to the 

• four-dimensional integral equations describing a system of three distinguishable 

relativistic particles. In this way we obtain the expressions needed to calculate 
all possible electromagnetic processes of the three-body system. An interesting 
^ , aspect of our results is the natural appearance of a subtraction term needed to 

^ ' avoid the overcounting of diagrams. 

> : 

•'-j . I. INTRODUCTION 

rS 

H , 

5t ! With the advent of new experimental facilities like the Thomas Jefferson National Accel- 

erator Facility (TJAF), the Electron Stretcher Accelerator (ELSA), and the Mainz Microtron 
(MUMI), there is currently great interest in the use of photons and electrons to probe the struc- 
ture of hadronic systems. In practice, this means using the electromagnetic probes to induce 
a variety of reactions amongst the hadrons. Here we shall be concerned with those reactions 
where effectively only one external photon is involved. This includes not only photoabsorption 
and photoproduction reactions, but also electron scattering and electroproduction when calcu- 
lated to lowest order in the electromagnetic interaction. On the theoretical side such reactions 
are described by {J'^), the matrix element of the electromagnetic current operator, and it must 
ultimately be the goal of the theorist to construct models where ( J^) describes the experimental 
data as accurately as possible. 

An essential constraint on (J^) is that it must obey current conservation, expressed by the 
continuity equation d^{J^) = 0. Current conservation is a consequence of charge conservation 
and is therefore a fundamental property of any theory. On the other hand, because one uses 
models to approximate the full theory, current conservation is not guaranteed a priori. For this 
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reason, much effort has been devoted to the question of how to impose current conservation 
within a particular model of strong interactions. In a seminal paper, Gross and Riska (GR) [|T| 
have shown how to construct the conserved current of two nucleons described by the Bethe- 
Salpeter (BS) equation. By analysing specific meson-exchange diagrams of the BS kernel, they 
showed that current conservation is achieved when the two-nucleon interaction current is con- 
structed by attaching photons to all possible places in the BS kernel. Although indispensable 
for the case of a two-body system, this result does not apply to systems consisting of three or 
more particles. There is also no straightforward way to use this result to construct the electro- 
magnetic current of three particles even in the case where only two-body strong interactions 
are present. 

Here we would therefore like to present a different method for constructing conserved cur- 
rents that is applicable to any number of particles whose strong interaction processes are 
described by relativistic four-dimensional integral equations. This method involves a direct 
gauging of the equations themselves in the sense that a vector index ^ is added to all terms of 
the equations in such a way that a linear equation in //-labelled quantities results. Proceeding 
in this way, we obtain integral equations for the gauged quantities of interest (e.g. the gauged 
Green function or the gauged scattering t-matrix T^) expressed in terms of other gauged 
quantities that are known or that can be easily constructed (typically the gauged one-particle 
propagator and the gauged two-particle potential f^). This approach, which we shall sim- 
ply refer to as gauging the equations, results in the external photon being effectively coupled 
everywhere in the strong interaction model, so that current conservation is guaranteed. 

Using this method for the case of two nucleons, we obtain that the hadronic current is a 
sum of matrix elements of the gauged nucleon propagators and the gauged BS kernel. Thus our 
result, in this case, coincides with the one of GR. However, our method is also easily applied to 
other systems. As the simplest strongly interacting system going beyond the two-body problem 
is that of three particles, we choose to illustrate the general nature of our gauging method by 
applying it to the relativistic three-body problem whose strong interactions are described by 
standard four- dimensional integral equations. That is, we consider systems like that of three 
nucleons or three quarks whose strong interaction processes do not involve coupling to the two- 
body sector. To keep the presentation as simple as possible we restrict the discussion to the case 
of three distinguishable particles. The gauging of three identical particles involves additional 
considerations of particle exchange symmetry and is the subject of the accompanying paper 0. 

Although four- dimensional three-body equations have already been pursued numerically by 
Rupp and Tjon [Q, there has not been a generalisation of the Gross- Riska result to the case of 
three particles. As a result, there is presently no rigorous derivation of the conserved current 
for a relativistic three-body systems. This paper is therefore devoted to showing how the 
gauging of relativistic three-body equations leads to gauge invariant expressions for the various 
electromagnetic transition currents of a three-particle system. It is a feature of our approach 
that the gauge invariance of our expressions is not imposed in an ad hoc fashion, but rather 
according to the way prescribed by field theory, namely by coupling photons to all possible 
places in the strong interaction model. 

To show the flexibility of the gauging method, we apply it to two different relativistic three- 
body equations. The first of these is the integral equation for the 3 — 3 Green function whose 
kernel is defined in terms of two-body potentials. This leads to three-body electromagnetic 
currents which are expressed in terms of two-body potentials v and gauged two-body potentials 
v'^. In the second approach, we gauge the four-dimensional version of the Alt-Grassberger- 
Sandhas equations 0] in order to obtain electromagnetic currents that are expressed in terms 
of two-body t-matrices t and gauged two-body t- matrices t^. Our final results consist of gauge- 
invariant expressions for the electromagnetic current of the three-particle bound state, as well 



FIG. 1. Examples of last cuts (dashed lines) in relativistic three- and four-body processes, (a) 
Ambiguous last cuts in the four-body problem, (b) Unique last cut in the three-body problem, (c) 
Ambiguous last cuts in electron scattering (top line) off three particles (bottom three straight lines), 
(d) The traditional picture of elastic electron scattering off a three-body bound state in impulse 
approximation. Interpreted as a Feynman graph (with the ellipses representing the bound state wave 
function), this diagram is consistent with the graph of (c) being included twice, once for each last cut 
shown. The diagram in (d) therefore contains overcounting. 

as the various electromagnetic transition currents of three particles. Thus, in the case of three 
distinguishable nucleons, our expressions describe the electromagnetic form factors of the bound 
NNN system (^H), and the processes Nd Nd'j, Nd 'jNNN, j^R^ Nd, and 7^11^ NNN. 

A notorious problem plaguing many four-dimensional approaches is the overcounting of di- 
agrams. For perturbation graphs, such overcounting can be corrected simply by subtracting 
the overcounted terms. However, when overcounting occurs within the scattering integral equa- 
tions themselves, the way to solve the overcounting problem becomes highly non-trivial. Indeed, 
only recently has such an overcounting problem been solved in the context of four-dimensional 
integral equations for the NN— nNN system P,^]. 

Overcounting can arise when reducible diagrams have ambiguous last cuts • An exam- 

ple of an ambiguous last cut for the four-body problem is given in Fig. |l](a). Luckily, there is no 
overcounting in the four-dimensional scattering equations considered here. This is because of 
the purely three-body nature of our system (there is no coupling to two-body channels) which 
makes all last cuts unique - see Fig. |l](b). However, as soon as coupling to an external photon 
is made, the system effectively obtains coupling to the four-body sector and overcounting again 
becomes a possibility. To see this explicitly, consider the calculation of electron scattering off a 
three-body system in the relativistic impulse approximation. One contribution to this process 
is given in Fig. 0(c). The diagram shown contains ambiguous last cuts of the same type as 
in Fig. |l|(a). This ambiguity needs to be carefully taken into account when summing the full 
perturbation series for the electron scattering process. For example, to find the expression 
describing electron scattering off a bound three-body system in the relativistic impulse approx- 
imation, one cannot simply use the direct generalisation of the non-relativistic expression as 
given diagrammatically in Fig. [l|(d). Interpreted as a Feynman graph, this diagram overcounts 
interactions between the two spectator particles (the lowest two lines in the diagram) since 
such interactions are contained in both the initial and final three-body bound-state wave func- 
tions. This overcounting corresponds to the diagram of Fig. |l](c) being included twice, once for 
each last cut shown. Just this type of overcounting appears to be present in Rupp and Tjon's 
calculation of the electromagnetic form factors of ^H. 

An important feature of the gauging of equations method is that it not only attaches photons 
everywhere in the three-body amplitude, but it also does this without introducing any over- 
counting. Indeed, it is found that the gauging procedure itself gives rise to subtraction terms 
that effectively remove all overcounted contributions. In this way the complications brought 
about by ambiguous last cuts like that of Fig. |l](c) are taken care of automatically. 

Preliminary results of the present paper (I) and the following paper (II) were first reported 
a few years ago 0. In the meantime we have appUed the gauging of equations method to 
the spectator formalism of Gross [|T^] to generate gauge invariant three-dimensional expressions 



describing any hadronic or quark system interacting with an external electromagnetic probe 



TT| . In particular, we have derived the three-dimensional expressions for the various electro- 
magnetic transitions currents of the three- nucleon system within the spectator approach [|l^. 
In this sense, the results of Ref. can be considered as a gauge-invariant three-dimensional 
reduction of the four-dimensional results presented in I and II. We have also applied the gaug- 
ing of equations method to the NN — nNN system where, as previously mentioned, the 
overcounting of diagrams provides an extra degree of complexity. For these works the present 
paper together with II form the basic theoretical foundation, and provide the references where 
all the missing details are given. 

Although the power of the gauging of equations method is well demonstrated on the example 
of the three-body system considered here, it should be emphasized that the same method is 
just as easily applied to other strongly interacting systems (quark or hadron) including those 
where the number of particles is not conserved. Indeed this method has recently been used to 
gauge the vrA^ system [14|, and as mentioned above, the NN— nNN system [p^] . 

As the gauging of equations method couples one external photon everywhere in the strong 
interaction model, it also forms the basis for exact descriptions of more complicated electromag- 
netic processes. For example, we have recently shown how to use this method to incorporate 



an mterna/ photon into all possible places in a strongly interacting system |T^. The resulting 
expressions provide a way to calculate the complete set of lowest order electromagnetic correc- 
tions to any strong interaction model described by integral equations. Being complete, these 
electromagnetic corrections are therefore gauge invariant. The gauging method can also be used 
to describe processes with more than one photon. For example, by gauging a strong interaction 
scattering equation twice, we would obtain gauge invariant expressions for the corresponding 
Compton scattering process. 

Finally, it is important to note that although we are concerned in this paper with the 
electromagnetic interaction for which gauge invariance (or current conservation) is a major 
issue, the gauging of equations method itself is totally independent of the type of external 
field involved. Thus all the expressions for transition currents developed in this paper hold 
also for cases where the external field is due to strongly or weakly interacting probes for which 
current is not conserved (of course the exact form of the gauged inputs would need to be chosen 
appropriately) . 



II. GAUGING TWO DISTINGUISHABLE PARTICLES 

A. Gauging the two-particle Green function G 

Gross and Riska |l| have shown how to construct a conserved current for the two-particle 
system described by the Bethe-Salpeter (BS) equation. In addition to the one-body current, 
which in general is not conserved, one also needs the two-body interaction current obtained by 
attaching a photon to each place inside all the Feynman diagrams defining the BS kernel. Such 
a current satisfies the two-body Ward-Takahashi (WT) identity [0 corresponding to the given 
model of strong interactions. On mass shell, this identity guarantees that the matrix element 
of the current operator satisfies the continuity equation. To obtain their result. Gross and 



Riska applied the one-body WT identities IjT^ to a number of meson exchange diagrams of the 
BS kernel. Here we would like to introduce a different method, that of gauging equations, by 
rederiving the Gross-Riska result for the case of two particles. In the subsequent sections, we 
shall illustrate the general nature of our method by applying it to the case of the three-particle 
system. 
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FIG. 2. Examples of the two types of Feynman diagram making up C^. (a) A diagram that can be 
constructed by attaching a photon to a strong interaction diagram of G. (b) A diagram that cannot 
be constructed by attaching a photon to a diagram of G. 

The two-body WT identity ||16| can be written as 

q^,G^\kik2]PiP2) = i[eiG{ki - q,k2;piP2) + e2G{ki,k2 - qipm) 

-G{kik2]pi + g,P2)ei - G{kik2;pi,p2 + q)e2] (1) 

where pip2 {kik2) are initial (final) momenta of the particles, the photon is taken to be incoming 
with momentum q so that ki + k2 = Pi+P2 + q, G is the full two-body Green function given by 

(27r)^5^(p'i +P2-P1 - P2)G{p[p2;piP2) = j d'^yid^y2d^xid'^X2 
and G^ is the corresponding five-point function 

^i{k,y,+k2y2-Pix,-p2X2) ({]\Tm'^^\y^)m^'^){y2)m'^^\xi)¥'^\x2)J^{'^)\Q). (3) 

In the above, and are Heisenberg fields of particle z, T is the time ordering operator, 
is the quantised electromagnetic current operator, |0) the physical vacuum, and Cj is the 

charge of the z-th particle. If the particles are isotopic doublets then Cj includes an isospin 

factor, e.g. for nucleons Cj = |[1 + T^\ep where is the Pauli matrix for the third component 

of isospin, and is the charge of the proton. 

For local field theory with distinguishable particles the electromagnetic current operator is 

given by 

n^) = E ^^^'^ + -^oVi^) (4) 

where ./other (^) consists of partial derivatives with respect to all other fields that are present in 
the Lagrangian C Note that we have used translational invariance to write Eq. (^) in terms 
of J'^(O), and we have defined both G and G^ to be without total 4-momentum conserving 5- 
functions. In Eq. ([^) we have used a notation where the two-body Green function G is labelled 
by four momentum variables even though only three are independent. Such notation allows 
us to write down the WT identity in an especially simple form that easily generalises to any 
numbers of particles. 

One can classify the contributions to G^' into two groups [|T^ . The first of these consists of 
all Feynman diagrams that can be constructed from G by attaching a photon to an appropriate 
diagram of G. An example is given in Fig. 0(a). The second group consists of diagrams that 
cannot be obtained from G by attaching a photon, an example of which is given in Fig. 0(b). 



The special feature of the second group is that each contributing Feynman diagrams satisfies 
gauge invariance all on its own. 

The goal of this paper is to show how to construct a gauge invariant by attaching a 
photon to all possible places in every Feynman diagram contributing to G. We shall refer to 
this attaching of photons everywhere as the gauging of G. On the other hand, diagrams like 
that of Fig. 0(b) are current conserving from the outset, and their construction is a separate 
problem which will not be considered here. Thus for the purposes of this paper, we take G^ to 
be the result of the gauging of G, with the understanding that the neglected contributions can 
always be added separately without affecting the question of gauge invariance. 

In a similar way, we use a superscript ^ to indicate the vector quantity obtained by at- 
taching photons everywhere to an amplitude or potential. However in this case we will require 
that no photons he attached to external lines of Feynman diagrams making up the amplitude; 
nevertheless, we shall still refer to the process of attaching the photons to all other places in 
the amplitude as gauging. 

The WT identity is usually given in terms of the exact strong interaction Green functions 
of the underlying field theory. It follows that this identity must also be valid to any order 
with respect the strong interaction coupling constant. Here we would like to point out that the 
WT identity is also valid for any particular single diagram of the strong interaction. This is 
because one can always construct a Lagrangian with respect to which the diagram in question 
represents the only case of some given order of the strong interaction. In this case, the G in 
the WT identity is given by just one diagram of the strong interaction, and G^ is the sum of 
the diagrams obtained from G by attaching the photon everywhere. In the same way, the WT 
identity is valid for any sum of diagrams included in G. Having this in mind, we consider G to 
be given by the Bethe-Salpeter equation where the kernel V is given by a model consisting of 
any number of connected two-particle irreducible Feynman diagrams. 

We start by expressing G in terms of its fully disconnected part Go, and the kernel V: 

G = Go + GqVG. (5) 

This is a symbolic equation that, for the case of two-particle scattering, represents a shorthand 
notation for 

(27r)^ (27r)4 ^°*^^'i^2' ^i^^Winr^, SiS2)G{siS2\ pm) (6) 

where it is understood that p'l + p'2 = Pi + P2 = + r2 = si + 82- In Eq. (P), neither the Green 
functions nor the kernel contain (5-functions corresponding to total momentum conservation. 
Thus the disconnected Green function Go contains only one 5-function and can be written as 

GoW2,PiP2) = {2nY5\p\-p{)dM)d2{P2) = {2n)'6\p', - p2)dM)d2{P2) (7) 

where di is the dressed propagator of particle i. To save on notation we write Eq. (0) symboli- 
cally as 

Go = d,d2 (8) 

where momentum labels and the momentum conserving 5-function [together with factor (27r)^] 
have been suppressed. 

Eq. (|^) is basically a topological statement regarding the two-particle irreducible structure 
of Feynman diagrams belonging to G; as such, it can be utilised directly to express the structure 
of the same Feynman diagrams, but with photons attached everywhere. Thus from Eq. (|^) it 
immediately follows that 



G'^ = + G^yc + Go^^G + CqVC^". (9) 

This result expresses in terms of an integral equation, and illustrates what we mean by 
gauging an equation, in this case the gauging of Eq. (||). Implied in Eq. (H) is the result 

[GaVGf = G^VG + GoVG + GqVG", (10) 

which illustrates a rule for the gauging of products that is rather reminiscent of the product rule 
for derivatives. Although Eq. (0) follows from a topological argument, we shall postulate such 
a product rule to hold also in other cases where topological arguments may not be applicable. 
Again, Eq. (||) is a symbolic equation representing 

(27ry^ (27ry^ ^i'^2)^('"i^2; SiS2)G{siS2; P1P2) 

+ J ^^^^Go{kik2;tit2)V''{tit2;uiU2)G{uiU2;piP2) 

(27ry4 (27r)^ ^0(^1^2; ViV2)V{ViV2]WiW2)G^{WiW2]PlP2) (11) 

where now the presence or absence of a photon with momentum q needs to be taken into 
account in specifying the momentum conservation relations: ki + k2 = Pi + P2 + q, where 
h + k2 = h + t2 = Vi + V2 = Wi + W2, and pi + P2 = ri + r2 = Si + S2 = Ui + U2- 

In Eq. (^, both G^ and Gq are obtained from the Green functions G and Go, respectively, 
by attaching photons everywhere. It is therefore important to note that the gauged potential 
V'^ is similarly obtained from V, but with no photons attached to external legs. This is because 
such contributions are already taken into account in the terms GqVG and GoVG^^ of Eq. (|[). 
Eqs. (H) and (|^) are a set of linear integral equations for G and G^, and could be solved as 
such if V and were given. However, we can also formally solve the equation for G^, and 
thus express it directly in terms of G. Simple algebra gives: 

G" - GoVG^" = G^ + G^VG + GoVG, 

G^ = {I- GqV)-^ [G^{1 + VG) + GoVG] = G (Gq ^G^Gq ^ + V^) G (12) 
where Eq. (|^) was used in the last step. Defining the electromagnetic vertex function by 

Eq. (1^) gives the essential result of this section: 

r'^ = r(; + (i4) 

where 



^0 — G'o ^GqGq \ (15) 



By gauging Eq. (H) one obtains 



G^ = dtd2 + (iirf^ (16) 



so that 



FIG. 3. Illustration of Eq. (14) expressing the two-particle electromagnetic vertex function as 
a sum one-body currents and the two-particle interaction current. 



r^ = r^d2i + drM (17) 

where is the one-particle electromagnetic vertex function defined by the equation 

df = d,r';d,. (18) 

Note that for nucleons, to lowest order in the strong interaction, rf = 6^7^ where 7^^ is a Dirac 
matrix. In Eq. ([T^ ) Fq is thus the sum of single particle currents, and is the two-body 
interaction current of the given model. Eq. (|14D is illustrated in Fig. |^. This is the essential 
result of Gross and Riska, and has been derived here using the gauging of equations method. 

B. Ward-Takahashi identity for 

The validity of Eq. (P) for the gauging of G is clear from the topological argument given 
above. We would nevertheless like to show explicitly that the C^, constructed in this way, 
satisfies the WT identity. To do this, we first prove that the gauged potential satisfies 
the WT identity even though no photons are attached to its external legs. Using a shorthand 
notation defined by 

we can write the WT identities for Gq and the quantity [GqVGoY 

2 

q^,G'^ikik2;piP2) = i^[eiGo{ki - q;piP2) - Go{kik2;Pi + q)ei\, (21) 

i=l 
2 

q^,[GoVGo]f'{k^k2;piP2) = tY,[(^iGoVGo{ki - q;piP2) - GoVGo{kik2;pi + q)e,]. (22) 

i=l 

Note that [GqVGo]'^ satisfies the WT identity because, unlike V^, it has photons attached 
everywhere. Although Eq. (|21|) and Eq. (|2^) would be automatically true if field theory were 
being solved exactly, here we work within a strong interaction model specified by the input 
quantities V and Gq] as such, these equations are assumed to be true by construction. 
Using the product rule of Eq. (p!oD , we have that 

q^.[GoVGo]''{klk2■,PlP2) = gMG'o^G'o(A;iA;2; piPa) + Co g^y^G'o(A;ifc2; P1P2) 

+GoVq^G^,{k,k2;pm)- (23) 



Writing out the integrals implied by this expression [see Eqs. (||) and (|TT])] and using Eqs. 
( P^ ) and (|^, we obtain that 

2 

q^y''ikik2]PiP2) = i ^ [eiV{ki - q; P1P2) - F(A;iA;2; Pi + q)ei] (24) 

4 = 1 

which is the WT identity for V^. We may now use Eqs. ( p4D and ( pi]) together with Eq. ([T^ ) 
to evaluate q^T^{kik2]PiP2)- Taking into account that — y = G^^, we obtain 

2 

q^jT^{kik2\piP2) = [G^'^ {kik2; Pi + q)ei - eiG'^ih - q;piP2)] ■ (25) 
1=1 

Using this result and Eq. (|13|) it immediately follows that 

2 

q^G^{kik2] P1P2) = i ^ [eiG{ki - q; P1P2) - G{kik2] Pi + q)ei] , (26) 
1=1 

thus proving the WT identity for the G^ obtained by the gauging of equations method. 

C. Gauging the two-body bound-state wave function ijj 

So far we have defined "gauging" to be the process where photons are attached to all places in 
perturbation diagrams. As Green functions and potentials have a diagrammatic interpretation, 
the gauging of these quantities has therefore a clear meaning. On the other hand, the bound- 
state wave function is a purely nonperturbative quantity, and thus cannot be associated with 
perturbation diagrams. One can nevertheless define the gauged bound-state wave function ip^^ 
by formally gauging the bound-state Bethe-Salpeter equation using our product rule. Thus by 
gauging the equation 

^ = G.Vi, (27) 

we obtain 

tl^^ = {GqVY ^ + G^ij^^ (28) 

in this way defining ip^. Here the symbol ip represents the two-body bound-state wave function 
ipp{piP2) defined by 

(27r)V(P-pi -P2)^p(j5ij52) = j d^xid^X2e'^P'^'+P''"'\0\T¥'\xi)¥^\x2)\P) (29) 

where = rn? with m being the mass of the bound two-body system. To save on notation, 
we do not write explicit spin and isospin labels in the state \P) nor in the wave function ipp^ 
nevertheless, such labels (and if necessary, sums over such labels) are to be understood as 
implicitly present. 



It is now straightforward to show that the '0'^ defined by Eq. (|28D is just the quantity 
ijj^{kik2] P) given by 



^^{kik2;P)= / (i^xirf^X2e'('=^"i+'=^"2)(0|r^(^)(xi)^(2)(^2)jM(0)|p) (30) 



where ki + k2 = P + q. The proof of this proceeds as follows. Assuming the field theory 
under consideration admits a two-body bound state, the Green function G, defined by Eq. 
exhibits a pole at = where P is the total four-momentum of the bound two-body system. 
Indeed from Eq. (|^) it can be shown that as P^ — > 

Gip,P2,PiP2) ~ t p2_^2 (31) 



where P = pi + P2 = p'l +P2 and -^p is given by Eq. p9|). Similarly from Eq. one can show 
that for — > m^, 

r^^^n u \ . ^^^(fcifcs; P)i)p{piP2) ,„^x 

G^[kik2,piP2) ~ t 55 5 • (32) 

with 'i/j^{kik2] P) given by Eq. (|5(]|). Using the last two results in Eq. (^), equating the residues 
at P^ = rn^, and writing ijj'^ as shorthand for ilj'^{kik2; P), we obtain that 

^A*^ = {GoVf + GoV^lj^i). (33) 

We therefore deduce that the given by Eq. (^) does indeed satisfy Eq. (pH]). Eqs. ( pTj ) and 
(|28| ) form a coupled set of equations which can be formally solved for tp^ giving 

r = {l~ G,V)-' {GoVY^ = GG,' {GoVfij 

= G{G,'G^,G,' + V^)ij. (34) 

Recalling Eq. (0) we obtain 

= GT^'ij. (35) 



Comparison with Eq. (|13|) shows that can be obtained from G^ by taking the right-hand 
residue at the two-body bound-state pole [this of course is also obvious from Eq. (0)]. As such, 
Gq^iP^ is just the transition current describing the photo/electrodisintegration of the two-body 
bound state, as discussed shortly. 

It is often convenient to work with the bound-state vertex function (j) defined by the equations 

^ = ^00 ; ^^ = 0^0. (36) 

To gauge <p we follow the same idea as for the bound-state wave function, namely, we define 0^ 
by formally gauging the bound-state equation for [which follows from Eq. (|27|) ]. In this way 
we obtain the coupled set of equations 

ct> = VGo(i) (37) 
r = {VGoY(l) + VGQr (38) 

which can be solved as above to give 

r = GtG{VG,Ycp. (39) 

Writing {VGoY = VGq + VG^ and G^^G = 1 + TGq where T is the t-matrix defined by 

G = Go + GoTGo, (40) 

Eq. (^) reduces after some simple algebra to the result 

0'^ = (i + TGo)r^V-r^V'- (41) 



D. Ward-Takahashi identities for and 0'^ 

The WT identity for ifj^ follows from Eq. (|35D and the WT identity for V^, Eq. (]25|). Writing 
Eq. (^) in its full numerical form 

r{hh;P) = J ^^^^G{hh;hh)T'^{hh;pmmPiP2;P)^ (42) 
we may use Eq. (|^) to obtain 

2 

^Yl [G'\ht2]Pi + q)ei - e,G-\ti - q-pip2)] ^p{piP2;P). (43) 



Since G ^i/jp = 0, the second term in the square brackets does not contribute and the above 
equation reduces down to the WT identity for i/j^, 

q^Vikih: P) = I [e^ki - Q, k^, P) + e^i^ik^, k^ - q; P)] . (44) 
The WT identity for 0^ is found in a similar manner. Written out numerically, Eq. ( P^f ) is 

(P^{k,k2;P) = J ^^[G'o^G](fciA;2;tit2)[^G'o]'^(tit2;PiP2)0(piP2;P). (45) 

Note that the square brackets here indicate the resultant function to which the following mo- 
mentum variables refer. As [FGo]'^ is a gauged input, it satisfies the Ward-Takahashi identity 
from the outset. We can therefore write 

q,r{k,k2;P) = J ^2n2r [Go'G]ik^h-,hh) 



2 

i=l 



{e^[VGo]{U - q;PiP2) - [VGo]{ht2■,p^ + q)e,} 0(pip2; P). (46) 



The first term in the curly brackets can be integrated over pi using the fact that VGo(f) = (p. 
Similarly the second term in the curly brackets can be integrated over ti using the fact that 
Gq^GVGq = Gq^G — 1. Only one term survives, 



q,r{k,k2.P) = I ^^XJ(27r)^5^(fc.-p.-g)e,</>(piP2;P), (47) 

which gives the Ward-Takahashi identity for 0^: 

q^rikik2; P) = I [e,<Pik, - q, k2; P) + e2<p{ki, k2 - q; P)] . (48) 

Note that neither of the WT identities of Eq. (M) nor Eq. ( ^S]) are zero, even for on mass shell 
ki and ^2. Thus neither nor (p^ obey current conservation. 



E. Two-body electromagnetic transition currents 



To obtain the physical t-matrix for any reaction involving an electromagnetic probe inter- 
acting with a quark or hadronic system, it is sufficient to specify the corresponding matrix 
element of the electromagnetic current operator. We refer to such a matrix element simply 
as an electromagnetic transition current and denote it using the symbol j'^. Thus, for exam- 
ple, a process like two-nucleon Bremsstrahlung NN 7AW is described by the AW — NN 
electromagnetic transition current Jqq given by 

J^o = G,'G^G,\ (49) 

The physical ^-matrix is then found by contracting Jqq with the photon polarisation vector e^. 

For processes with one two-body bound state, like deuteron photodisintegration, one de- 
termines the electromagnetic transition current Jq from by taking the appropriate residue 
at the bound-state pole. Thus, to determine the {ij) —>■ ij electromagnetic transition current 
describing the process j{ij) ij where a photon is absorbed on the bound state (ij) producing 
free particles i and j in the final state, we take the residue of Eq. (0) on the right to obtain 

jI; = G.'GT^ij (50) 
= (l + TGo)r^V' (51) 

where T is the off-shell two-body t-matrix. As noted previously, a comparison with Eq. (H) 



shows that the (ij) ij electromagnetic transition current and the gauged bound-state wave 
function are simply related by 

jI; = G,'r. (52) 

It is also interesting to examine the relationship between j^ and the gauged bound-state 
vertex function 0^. From Eq. (^) one immediately obtains that 

j^ = <P^ + T^,Go<p. (53) 

From this equation it is evident that the difference between 0^ and the electromagnetic transi- 
tion current j^ is that 0^ has no photons attached to external legs. This suggests that <p'^ does 
not conserve current, while j^ does. Indeed q^j^ = G^^Qfj^ip^ = for on mass shell external legs 
since the momentum shifts {—q) contained in the WT identity for ip'^, Eq. (0), ensure that the 
poles of qfj,ip'^ do not cancel the zero of Gq^. 

Lastly we consider the case where an electromagnetic probe interacts with two particles 
which are bound in both initial and final states, for example elastic electron-deuteron scattering. 
Such processes are described by the electromagnetic bound-state current j'^ which can be found 
from Eq. (|T3|) by taking residues at both the initial and final bound-state poles. As these bound 
states can have different total momenta, we write the bound-state current as 

f = ^kT^^p (54) 

where = = m^. Current conservation for j^^ is obtained by using the WT identity for 
r/^, Eq. (^, and noting that ^G-^ = G-^ijj = 0. 



III. GAUGING THREE DISTINGUISHABLE PARTICLES 



Having gauged the two-body system, we are ready to apply the gauging of equations method 
to the substantially more complicated case of three strongly interacting particles. 

A. Gauging the three- particle Green function 

The Ward-Takahashi identity for three particles can be derived by following the same pro- 
cedure as used by Bentz [TO for the two-particle case. One obtains 



q^,G''{kik2h;piP2P3) = i[eiG{ki - g, /i;2fc3; + &2G{kx, ^2 - g, k'i;pxV2Vi) 
+e3G{kik2,k3 - q]pip2Pz) - G{kik2k3]pi + q,P2Pz)ei 
-G{kik2k3;pi,p2 + q,P3)e2 - G{kik2k3; pip2,P3 + q)e3]. (55) 

By extending the shorthand notation of Eqs. ([19[) and (|20|) to the case of three particles, the 
three-particle Ward-Takahashi identity can then be written more concisely as 

3 

q^,G^{kik2k3;piP2P3) = i^[eiG{ki - q;piP2P3) - G{kik2k3;pi + q)ei\. (56) 

1=1 

Here G is the full Green function of the three-particle system 

{27Ty6^{p[ + P2 +P3-P1-P2- p-i)G{p[p'2P3: P1P2P3) = j d'^yid'^y2d^y3d^xi(fx2d'^X3 

e*{pi2/iMy2+pJ,y3-Pixi-p2X2-p3X3) (l^\T^(^\y{)^^'^\y2)^'^^\y3)^^^\x{)¥'^\x2)^^^\x3)\Q) 

(57) 

and G^ is the corresponding seven-point function 

G^{kik2k3;pip2P3) = J d'^yid'^y2d'^y3d'^Xid'^X2d'^X3 

-p-^--^){0\T¥^\yi)¥^\y2)¥^\y3)¥^\xi)¥^\x2)¥^\x3)J^{0)\0). 



J{kiyi+k2y2+k3y3-piXi-p2X2- 



(58) 



In practise, the Green function for three distinguishable particles is specified by the integral 
equation 

G = Go + GqVG (59) 

where the potential V is three-particle irreducible, and in the absence of three-body forces, is 
written as a sum of three disconnected potentials V^: 

V = Vi + V2 + V3. (60) 

Here we use the usual spectator notation of three-body theory: defining {ijk) to be a cyclic 
permutation of (123), Vi is the potential where particles j and k are interacting while particle 
i is a spectator. Explicitly we have that 

yi{p'iP2P'3^PiP2P3) = Viip'jp'k, PjPk)d~^ (pi) (27r)^5^(p. - Pi) (61) 



where f , is the two-body potential between particles j and k. Neither Vi nor Vi contain total 
momentum conserving ^-functions, and the presence of 6^{p[ — Pi) in Eq. (^) is due purely 
to the disconnectedness of the potential Vi. In keeping with our shorthand notation, we write 
Eq. (0) as 

Vi = Vidr\ (62) 

Our goal is to gauge Eq. ( |5PD and in this way obtain the electromagnetic transition currents 
for the three-body system. As the gauging procedure is symbolically identical to that already 
presented for the two-body system, we conclude that Eqs. (|TB|)-(|T5|) apply also to the three-body 
system. Thus the seven-point function can be written as 

= GT^G (63) 

where the electromagnetic vertex function for three particles is given by 

= Gq^G^Gq^ + v. (64) 

In the three-particle case, however, V is given by Eqs. (|60|) and ( |62D so that 

V^' = V^^ + V^ + V^. (65) 

To obtain Vf^ all that is needed is to gauge the Vi of Eq. (^). Yet because of the presence 
of the inverse propagator in ^, there is some question of how to do this consistently. An 
unambiguous answer is provided by gauging Eq. (|59D with V exphcitly given by Eq. (|6^) and 
Eq. (^). In this case we write Eq. (^) as 

G = Go + 5^ G^VG = Go + 5^ D^^ViG (66) 

i i 

where Doi is the fully disconnected two-particle propagator defined by 

Doi = djdk- (67) 

The gauging of Eq. (^) thus involves the term 

{GoVGr = {Do^v.Gr = {D^,v, + D,,vt) G + D,,v,G^. (68) 

Now using the fact that D^^ = (Gg — Doid^ )d^^, the above equation gives 

(GoVGr = [Gi;V - Corf t;, + Go {v^d;')] G + Go^.G^^ 

= G^oViG + Co (-rf + v^d;^) G + Go^.G^ (69) 

from which it follows that 

= v^d-' - v.vl (70) 

It is seen that the same result would follow from gauging Eq. (|62|) directly using our product 
rule, as long as we use the prescription 

K')" = -rr. (71) 

This prescription can also be obtained by formally gauging the identity did^^ = 1 and taking 
1^ = 0. The expression for the gauged potential Vf", Eq. ([fOD, is illustrated in Fig. ^. The 



FIG. 4. Illustration of Eq. (l70|) for the gauged disconnected potential V-^ . 



negative sign in front of the term ViV^ may appear to be surprising, yet it is just what is 
needed to stop overcounting. Consider for example the gauging of the term GqVG appearing 
in Eq. (||): {GoVGY = G^VG + GqV^'G + GqVG''. It is apparent that the rightmost diagram 
of Fig. H appears in each of the three terms GqVG, GqV^G, and GqVG^. Thus the negative 
sign in question is needed to ensure that this diagram contributes only once to the gauging of 
GoVG. 

The fully disconnected three-particle Green function Go is given by 

Go = didjdk (72) 

where two momentum conserving ^-functions are implied in the used shorthand notation. Gaug- 
ing this expressions gives 

3 3 

G^ = 5Z<^.4 = 5Z</^oi (73) 

i=l i=l 

where the sum is over the three cyclic permutations of (ijk). The impulse term in Eq. ([6^ ) 
can then be written as 

Go'Gi;G,' = f2^fD,^. (74) 

i=l 

Using this and Eq. ( |7DD in Eq. ( |51D then gives 

= E (rf^o.^ + v^dr' - v^rt) (75) 

i=l 

which we illustrate in Fig. p. 

Noting that the full two-body Green function Di satisfies the equation 

A = + D^iViDi, (76) 

we can also write Eq. (^5]) as 



r 



+ 



V 



V 



FIG. 5. Illustration of Eq. ( |75D expressing the three-particle electromagnetic vertex function F'^ 
as a sum of one- and two-particle currents. 



3 

= E (^'A^' + ■ (77) 

i=l 

Eq. (^) [or equivalently Eq. (^)] is the main result of this section and extends the Gross- 
Riska resuh of Eq. (|T^) to the three-particle sector. It gives the precise way that the one- and 
two-body currents need to combine in order to obtain proper gauge invariance. 



B. Three-body bound-state current 

The three-body bound-state current describes processes where a system of three strongly 
interacting particles is bound both before and after the interaction with an electromagnetic 



probe. Examples include elastic electron scattering from H and He |]T9[, and elastic electron 
scattering from the proton considered as a bound three-quark system |^ . 

The bound-state current is found by taking left and right residues of Eq. ( |63D at the 
initial and final three-particle bound-state poles. With M being the mass of the bound state, 
it can be shown that 

G{PiP2P^'^PiP2Pz) ^ ^ p2 _ ^2 P (78) 

where \E'p is the three-particle bound-state wave function defined by 
{2nY5\P - pi - P2 - p-i)"^ p{piP2P-i) = 

/d*...*..A3e-''-™>(0|«".,.0*'^'fe)*'»'(.3)|P). (79) 

Here \P) is the eigenstate of the full Hamiltonian corresponding to the three-particle bound 
state with momentum P^. We thus find that 

= qjj^T^'^p (80) 
= j2^K (rf^oY + - v.r^) ^P. (81) 

i=l 

Note that the wave function \E'p used here satisfies Eq. (^), and thus fulfils the normalisation 
condition (for distinguishable particles) 

'^P^2 (Go' -V)^p = 1. (82) 

This result follows upon taking residues of the identity G = GG~^G and using the fact that 
G~^ = Gq^ — V. By exposing the bound-state poles in the G^ of Eq. (pSj), one finds that is 
also the matrix element of the current operator 

f = {J^^) = {K\,nO)\P) (83) 

where = P^ = M^. In terms of j'^, the familiar charge and magnetic form factors of spin-1/2 
bound three-body systems are given by 

FUf) = E (8") 

(86) 

where ^ indicates a sum over final spins and an average over initial spins (recall that states 
\P) and wave functions are to be understood as having implicit spin and isospin labels). In 
Eq. (p5D n is the magnetic moment of the bound state, Mj^ is the nucleon mass, and = —q^. 



1. Role of subtraction term in the bound-state current 



When used to calculate the bound-state current j'^, the role of the subtraction term in 
Eq. ( [75|) is seen especially well by writing the bound-state wave function in terms of its usual 
Faddeev components: 

3 

^ = ^^„ (86) 

i=l 

where \l/ is written without a momentum subscript to save on notation, and where 

= Do^v,^. (87) 
For this purpose it is sufficient to consider just the one-body current contributions to j^, 

3 

irne-body = E'^^'(^o/-^0^- (88) 

1=1 

Using Eq. (^) this then evaluates to 

3 

Jo'ne-body = E ^^'D^'C^, + ^ ,) . (89) 

i=l 

This equation is not symmetrical with respect to the initial and final state wave functions. 
That this must be so can be seen as follows. The term r'^D^^ corresponds to a first order 
electromagnetic interaction with particle i. During this interaction, particles j and k can be 
interacting, and this contribution is contained fully in the final state wave function ^. To avoid 
overcounting, there must be no preceding j — k interaction coming from the initial state wave 
function. This indeed is the case since, as is evident from Eq. (p7|), the very last interaction in 

is between particles j and k, and just this component is missing from the right-hand side 
(RHS) of Eq. (|9|). 

Another way to write Eq. (|89D is 

3 

Jo"ne-body = ^^'o^^ " E ^o/^. (90) 

1=1 

where Fq = Yl^=i ^i^o-^ electromagnetic vertex function in the impulse approximation. 

In a three-dimensional approach, the term 'i'FQ\E' would correspond to the correct expression 
for the total one-body current. In a four-dimensional approach, however, one needs to subtract 
the last term of Eq. ( pPj ) in order to avoid the overcounting just discussed. In this regard we 
note that the pioneering four-dimensional calculation of Ref. used ^^Fq^ as the expression 
for the one-body current, and therefore the results of this calculation contain overcounting. 

Further insight into the nature of the overcounting problem can be gained by assuming, as 
was done in Ref. [Q, that the input two-body potentials are separable. In this case we have 

Vi = hiXhi (91) 

where hi is the separable potential form factor with particle i being a spectator. The two-body 
t-matrix is then given by 



ti = hiTihi (92) 

where 

r-' = A-^ - hiD^ihi. (93) 

For separable potentials we follow Ref. and introduce the "spectator bound-state form factor" 
Xi defined by 

= GohiTiXi. (94) 
As the Faddeev equation for the wave function components is 

= Do^U (^, + ^k) (95) 

where ijk are cyclic, it follows that 

Xi = dr'hi (^,- + ^fe) . (96) 
Using these relations in Eq. (R^) it is easy to show that 



one-body 



J2 i X^d^nX, + X,r,h,Do^dthkTkXk j (97) 

which has a straightforward graphical interpretation. Depicted in Fig. ^(a) are the topologically 
distinct contributions to the sums in Eq. (p7|). One can similarly write 



(9J 



which compared to the previous equation, contains an extra term given by 

^r>;V,^ = Y,X^dtT^K\^X,. (99) 

i 

This term, depicted in Fig. ^(b), represents what is being overcounted in Eq. (^). Considering 
this overcounting term together with the first term in Eq. (p8| ) [depicted by the first term 
in Fig. |^(a)], we see that the role of the overcounting term is to overdress the quasiparticle 
propagator while the spectator particle is being gauged. Thus, in the separable case, the use 
of the expression ^'rQ\l/ to calculated the one-body electromagnetic current is equivalent to 
overdressing the quasiparticle propagator in the first term of Fig. ^(a). 



(0) ^if^ + ^r^^ - 





(b) 

FIG. 6. (a) The one-body electromagnetic current when two-body interactions are separable, (b) 
The overcounting term present in ^'Fq'I'. The double- line with the black dot represents the dressed 
quasiparticle propagator r, the small open circle represents the separable potential form factor h, and 
the large ellipses represent the spectator bound-state form factor X. 



2. Gauge invariance of the bound-state current 

To prove the gauge invariance of Eq. (^) for the three-body bound-state current, it is 
convenient to use a symbohc notation for the WT identities. To do this we introduce the 
operators Cj whose numerical form is defined by 

ei{kik2h,piP2P3) = iei{2'Kf'^5'^{ki- Pi - q)b^{kj - Pj)S'^{kk - Pk) (100) 

where ijk represents a cychc ordering of 123. Then the WT identities for the gauged two-body 
potential and gauged one-particle propagator can be written in three-particle space in terms of 
commutators as 

q^.v'^Ii = [cj + 4, Vi], q,,d''jjlk = [e^, di] (101) 

where Jj, Ij, and are unit operators in the space of particles i, j, and k, respectively. 
Writing Eq. ( |75D in the form 

3 

= E {di'dfdr'dj'd,' + v^dr' - v^d-Xd^') (102) 

i=l 

and using the WT identities of Eq. ( |101|) , the first term on the RHS of Eq. ( |102|) gives 

3 

q, J2 di'd^d-^'d^' = -[e, G,'] (103) 

i=l 

where 

e = 61 + 62 + 63. (104) 

Similarly for the last two terms of Eq. ( p.02| ) we have that 

q^ {v'j^d''^ - Vidi^d'^d^^) = [ej + ik, + M^i^ d^^] = [e, Vid^'^] (105) 

and therefore 

3 

9m E K^r' - ^^dr'€d^') = h V] (106) 

i=l 

where V is the full potential as given by Eq. (|60D. We have thus shown that 

g^r^ = -[6,^-1] (107) 



where G'^ = Gq^-V is the inverse of the full Green function. It is recognised that Eq. ( [107| ) 
is just the shorthand three-particle version of the two-particle result given in Eq. (|25|) . Current 
conservation and therefore the gauge invariance of the bound-state current follows immediately. 
It is worth noting that the presence of the subtraction term [last term in Eq. (|102|) ] is essential 
for current conservation. 



3. Generalization to all transition currents 



Electromagnetic currents for all possible transitions can be constructed by following the 
same procedure as above for the three-body bound-state current. That is, one can begin with 
the expression = GT^G and then take residues at the relevant initial and final state poles. 
In this respect it is important to note that the Green function G has poles corresponding not 
only to single particle propagators and three-body bound states, but also to two-body bound 
states. For example, it can be shown that in the vicinity of a two-body bound state of particles 
2 and 3 

G{k,k,h,p,p,p,) ^ ^ ^^-^t'^^^^P^^ n"V(PiP2P3) as K'.^m' (108) 

where m is the bound-state mass of particles 2 and 3, and Ki = k2 + k^. Similarly 

G{kik2k,,pip2P3) ^ in,p,{hk2k3) '^^^]lf'^^^f'^ as P^^m' (109) 

where Pi = P2 + Ps- The seven-point function G^ has the same poles as G and one can show, 
for example, that in the vicinity of these poles 

r^un 1 1 \ .■ipKi{k2k3)d{ki) i iM/nM d • \ -^bO^A fePs) 

G^{kik2ks,pip2P3) ~« 5 out,iriA;i J''(0) Pipi,m)z — . (110) 

Kf — irt'^ Ff — 

Analogous expressions for G and G^ hold in the vicinity of any poles of G corresponding to 
physical initial and final states. The in- and out-wave functions appearing in these expressions 
are defined generally by 

{2nY5\P -PI-P2- P3)^f (P1P2P3) = 

j d'^xid^X2d^X3e'^P''''^P''''^P''''\0\T¥^\xi)¥^\x2)¥^\x3)\jAn), (111) 

i27Ty6\K -ki-k2- k3)^r\hk2k3) = 

d^yirf^y2A3e-'('^^^+'^^^+'^^^)(out,z|T^'(i)(yi)^'(2)(2/2)^^'H?/3)|0), (112) 



where |j, in) and |i,out) are the in- and out-states of quantum field theory (QFT) for the 
physical initial and final states of total momentum P and K, respectively. Here we follow the 
usual convention where index j = 1,2,3 labels a state where particle j is free and the other 
two particles are bound, and j = labels the state where all three particles are free. In case 
of an initial three-body bound state, the index j and label "in" are dropped (for three-body 
bound states there is no difference between in- and out-states as only one physical particle is 
involved). Similar conventions hold for index i. For a given model specified by the three-body 
potential V, the in- and out-wave functions are specified by the equations 



= ^f) + ^°^t(X/-\/^)G', "^f = '$f + Gj{V -Vj)¥l' (113) 

where for z = 0, = 0, "^f^ is the wave function of three free particles, and Gi is the free 
Green function, while for i = 1, 2, 3, is given by Eq. (|62D, '^f'^ is the wave function of a free 
particle i and a bound [jk) pair, and Gi is the disconnected Green function with particles j 
and k interacting and i being a spectator, i.e. 



Gi — diDi. 



(114) 



Eqs. ( |113|) hold also for three-body bound states if we drop the inhomogeneous terms '^f^ and 
'^''^^ but otherwise use i,j = 0. 

Taking right and left residues of at the initial and final state poles, it is easy to see from 
the above equations that the general expression for an electromagnetic transition current is 

j^;. = (out,^|J^(0)|j,in) = ^routp^vl/;-. (115) 

Eq. (|115|) is valid no matter whether the initial or final states consist of three free particles, 
one free and two bound particles, or the three-body bound state. Likewise the proof of current 
conservation for Eq. (|1 1 5|) is the same for all cases: one uses Eq. ( |107|) and the fact that 



= 0; G-^^f = (116) 



to deduce 



d'^ki d^k2 d^pi d^p2 



^ I (27r)4 (27r)4 (27r)4 (27r)4 

3 



iJ2^r\kik2h) [G-\k^k2k^,pi + q)ei - eiG-\h - q\ViV2Vz)\ ^^{vmVz) = 0. (117) 



1=1 



C. Gauging the AGS equations 

In the previous subsection we constructed the electromagnetic currents for all possible tran- 
sitions of the three-particle system. The expression derived, Eq. (|115| ), expresses the transition 



currents in terms of the vertex function F'^ which in turn is given in terms of the gauged 
three-body potential - see Eq. (p^). Although this formally solves the problem of how to 
gauge a three-body system, it needs to be recognised that Eq. ( |115| ) may not always be very 
useful for practical calculations. For example, in the present four-dimensional case where the 
three-particle potential is assumed to be of the simple form given by Eq. (|60D and Eq. (|62D, 
the numerical evaluation of Eqs. (|113|) for the scattering in- and out-wave functions is problem- 
atic. This difficulty is due to the disconnectedness of the three-particle potential which forms 



the kernel for these integral equations. For other cases the use of Eq. (|115|) can similarly be 
impractical. In the four-dimensional description of the nNN system the possibility of 

pion absorption together with overcounting problems makes the tiNN potential V difficult to 
specify, let alone calculate. Likewise in the spectator approach to the three-nucleon problem 
where spectator nucleons are put on mass shell , the effective three-nucleon potential is not 



easily revealed. 

All these problems can be resolved by avoiding the use of potentials in the formulation 
of the strong interaction three-body problem. This is just what is done in the case of three- 
dimensional Quantum Mechanics by using the Alt-Grassberger-Sandhas (AGS) equations to 
describe the strong interactions of three particles. The AGS equations take as input two-body 
t matrices, rather than potentials, result (after one iteration) in a connected kernel, and form 
one of the standard tools for doing practical three-body calculations. These same benefits can 
be obtained in four dimensions by using equations that are the four-dimensional analogue of 
the AGS equations. Indeed such equations (which we also refer to as the AGS equations) have 
already been used in the successful formulation of the nNN problem [§,0, and in the three- 
nucleon problem within the spectator approach [jlO[. Thus for practical reasons it is important 



to apply our gauging of equations method directly to the AGS equations themselves. In this way 
we shall obtain a gauge invariant generalisation of the AGS formulation to systems consisting 
of three particles with an added external photon. 



1. j{ki) i{jk) transition current 

Our starting point is the four-dimensional version of the AGS equations describing the 
scattering of three strongly interacting particles. These AGS equations can be written in the 
two forms 

3 3 

Uij = Gq ^bij + dikT^GoUkj ; Uij = Gq ^6ij + UikGoT^Skj (US) 

k=l k=l 

where the AGS amplitude Uij describes the process j{ki) — > i{jk), i.e. the scattering of particle 
j off the (ki) quasi-particle, leading to a final state consisting of particle i and the (jk) quasi- 
particle. If the quasi-particles form bound states then the amplitude Uy is related to the 
t-matrix for the physical j{ki) i{jk) process by 

T,, = i^iUjiPj (119) 

where i/ji is the two-body bound-state wave function of the (jk) system. In Eq. ( |118|) , is 
defined by 

n = hd^' (120) 

where t^. is the two-particle t-matrix for the scattering of particles i and j. The AGS equations 
of Eq. ( |118[ ) can be written in matrix form as 

U = Gq^I + irGoU ; U = Gq^I + UGoTI, (121) 

where the [i, A;)'th elements of matrices U, T, and X are defined by 

Mik = Uij, (122) 
[Uk = SikT,, (123) 

[^L = ^^k = l- (124) 

Although one could now gauge the matrix U by gauging Eqs. (|121|) in the usual way, the 
presence of Go ^ in Eqs. (|T^ ) makes it more convenient to instead gauge the Green function 
quantity 

U = GqUGq (125) 

which satisfies the equations 

U = IGo+IGoTl{ ; U = GoI + UTGoT. (126) 

At this stage it is convenient to make use of the two-body bound-state vertex function 0j, 
defined as in Eqs. (^) by 

iJi = Doi(f)i ; tf^i = (l)iDoi. (127) 
Then the t-matrix of Eq. ( |119[ ) can be written as 



Tij = (piDoiUijDoj(j)j = (j)idi ^Uijdj Vj- (128) 

As previously discussed, gauged potentials and t-matrices do not have photons attached to 
external legs. Yet such attachments are necessary for gauge invariance. For this reason we do 
not gauge the t-matrix Tij, but instead introduce the quantity 

fij = d^Tijdj (129) 

which contains extra propagators for the initial and final state spectator particles. It is by 
gauging Tij that photons get attached to all possible places in the j{ki) — > i{jk) process. 
Of course after gauging, it is necessary to remove these propagators to get the corresponding 
electromagnetic transition current jij, thus 

= d;'fi;dj\ (130) 

From Eq. (|T28|) we find that 

Tij = (f)iUij(pj. (131) 

Gauging this equation gives 

Ttj = '^tU^^J + ^^Ui,r, + ^^^4>y (132) 

Here 0f and 0^ are the gauged two-body bound-state vertex functions discussed in Sec. II. 
Being gauged quantities of the two-body problem, they form an input to the gauged three- 
body problem. As Uij is assumed to be known from the solution of the strong interaction 
three-body problem, only the gauged AGS Green function is left to be determined. We 
could find by gauging Eqs. ( |126| ) explicitly; however, this is not really necessary as there is 
a one-to-one correspondence with the previous gauging of Eq. (|59|) . Indeed Eq. ( p.26| ) follows 
formally from Eq. ( pD| ) upon the following substitutions: 

G^U- Go^IGo; V^T. (133) 

Moreover, just as Eq. (|60|) expresses as a sum of three components Vi = Vid^^, we can 
similarly write T as a sum 

r = Wi + W2 + W3 (134) 

where 

Wi = wdi^. (135) 

As matrix T is specified by Eqs. ( |12(]| ) and ( |123| ), it follows that Wi is a matrix whose (n, m)'th 
element is given by 

['f^i]nm — ^ni^imti- (136) 

With the correspondence now complete, we can immediately use Eq. ([75| ) to write down the 
gauged matrix lA^. We obtain 

W = UVU (137) 
where F'^ is a matrix whose (n,m) 'th element is 



3 

^'nm = E i^nL^iD,- + ^m^^mt^dT^ - J^.^.^t.rf) (138) 

4 = 1 

3 

= ^'^o' + ^n^mKd-' - 6nmU^. (139) 

1=1 

Alternatively, we may write in matrix form as 

= E i^~'^"Do' + - ^.rf ) . (140) 

The use of the same symbol for both the vertex function of Eq. ([75| ) and the matrix of 
Eq. (|140|) should not cause confusion as only the latter appears in matrix expressions. On 
the other hand, using the same symbol has the advantage of emphasising the formal similarity 
between the two. For example, the matrix form of Eq. ( |14(]| ) is equally well illustrated by Fig. ^ 
with v^s replaced by w's. 



2. Current conservation 



In this subsection we would like to show that the current j^j as specified by Eqs. ( |130D and 
( |132|) is conserved. To do this we first show that the gauged AGS Green function satisfies 
the Ward-Takahashi identity. Writing Eq. ( |140| ) as 

3 

= l-^G,^G^,G,^ + E K^r' - «^.rf ) , (141) 

i=l 

we may use the Ward-Takahashi identities for the input quantities ■wf and : 

q^w'^{kjkk;pjPk) = i'Y[eiWi{ki - q;pjPk) - Wi{kjkk;pi + q)ei] (142) 
q^v'^ih^pi) = i [d~^{ki)ei - eid:r^{pi)] , (143) 



where the former equation follows from identical arguments to that proving Eq. (|24f) , and where 
the latter equation follows from 

qf,d'^iki,pi) = i[eidi{pi) - di{ki)ei]. (144) 

From the Ward-Takahashi identity for Gq, it is also easy to show that 

3 

q^r'^{kik2k3;pip2P3) = [Gq ^(^i^a^s; Pi + q)ei - CiGg'^iki - q;piP2P3)] (145) 

i=l 

where 

r[; = G,'G^,Go\ (146) 
Using these results in Eq. ( |141|) gives 



qi,T''{kik2h;piP2P3) = {X ^[GQ^{kik2h;pi + q)ei - eiGQ^^h - q;piP2P3)] 



1=1 



+d. ^iki,pi) "^[eiWiik - q]PjPk) - Wi{kjkk\Pi + q)ei] 
-Wi{kjkk;PjPk)[di\ki)ei - eid^^Pi)]} 



(147) 



As ZlLi Ei^i (^ih = ELi Ez^i we can rewrite Eq. ( |Ii7| ) as 
qf,T^{kik2k3;pip2P3) = 

I'^GQ^{kik2k3;pi + q) -^dY^{ki,pi)wi{kmkn;Pi + q) - Wi{kjkk\PjPk)dj^{ki 



i=l 



i- 

3 



4=1 



J ^^^^(/Ci - q;PiP2P3) - X^c?; ^{kupi)wi{ki - q\PmPn) - Wi{kjkk; PjPk)di ^{pi 



(148) 



where iijk^ and {Iran) are both cychc permutations of (123). Recognising that the square 
bracket terms of Eq. ( |148 ) correspond exactly to the expression for lA^^ derived from Eqs. (126), 
we may write 



6 

qiiX^{kxk2k-i\pxP2P3) = * \U'^{kik2kj,\Pi + q)ei - tiU^^iki - q;piP2P3) 

i=l 

Using Eq. ( p.37] ), we thus obtain the Ward-Takahashi identity for W^: 



(149) 



q,M^{kik2k^]PiP2P3) \eiU{ki - q;piP2P3) -U{kik2k^]Pi + q)ei . (150) 

i=l 

Next we use this result and the last term of Eq. (|132|) to write 

qti(t>iUtj4>j{kiKi;pjPj) = i / {2TTy^5^{Ki - kj - kk)5'^{Pj - Pi - Pk)d'^kj d^kk d^pk d^pi 



i{kjkk] Ki) \^iUij{ki - q;piP2P3) - Uij{kik2k^]Pj + q)ej (pjipkPu Pj) 
+ ^ (piikjkk, Ki)eiUij{ki - q;piP2P3)<Pj{PkPu Pj) 



- ^ (t)i{kjkk\ Ki)Uij{kik2ks;pi + q)ei(j)j{pkPi; Pj 



(151) 



where the momentum variables are as specified in Fig. |^ To find the contribution of this term 
to q^Jij, we need to multiply it by d~^{ki)dj^{pj) and then take the on mass shell limit, i.e. 
d^^{ki) and d~^{pj) — > 0. Doing this we see that the first term in the curly brackets will 
give zero and in this sense is gauge invariant. On the other hand, the last two terms in the 
curly brackets will not give zero since the factor d^^{ki)d~^{j)j) will be cancelled by propagators 
di{ki) and dj{pj) contained in the Uij. Thus, the contribution of the last two terms of Eq. ( [151| ) 



to jl^j will not be gauge invariant. However, it is now easy to check that gauge invariance is 
restored by including the other two terms defining Tj^j in Eq. ( |132| ). To do this, we change 
integration variables in the last two terms of Eq. (|151|): 





FIG. 7. 
Eq. {^). 



Illustration of the variables involved in the evaluation of q^(j)iUj^j4>j{kiKi;pjPj) as in 



qf,(f)iU^j(l)j{kiKi;pjPj) = i / d'^kj d'^kk d'^pk d^Pi {2n) ^ {5'^{Ki - kj - kk)5'^{Pj - Pi - Pk) 



(piikjkk, Ki) eiJJij{ki - q; pmPs) - Uijikik2k3; pj + q)ej (pjiPkPf, Pj) 

+6'^{Ki - kj ~kk- q)5'^{Pj -Pi- Pk) 4>i{ki + q; Ki)eiUij{kik2k3;piP2P'i)(t)j{pkPu Pj) 

-5^{Ki - kj - kk)5^{Pj + q-Pi-Pk)^ Mkjkk, Ki)Uij{kik2ks;piP2P3)ei(j)j{pi - q; Pj 



(152) 



Making use of the WT identity of Eq. (0), we have that 



%(l>^j{VkVh Pj) = i ^ ei(f)j{pi - q; Pj 

1^3 



J/' 



q^,(j)fikjkk; Ki) = ^^^^^ + ^' ^^^'^^ 



and it becomes evident that the last two terms of Eq. (|152|) correspond exactly to —q^(f)^Uij(j)j 
4>iUijq^(pj. Thus contracting Eq. ( |132D with q^ gives 



q^Tt^ = I J d% d^kk d'^pk d'^p, {2n)-''S\K, - k, - kk)S\P, - - pk) 
(j)i{kjkk] Ki) eiUij{ki- q;piP2P3) -Uij{kik2k3;pj + q)ej (f)j{pkPi; Pj), (154) 
and the current conservation of r^- follows. 



3. Unwritten without subtraction terms 



The gauged AGS Green function [/^^ is expressed by Eqs. ( [137|) and ( |139| ) as 



i=l 



(155) 



This form for U^^^ may be the best for practical calculations, but the presence of the minus 
sign in the subtraction term —Uniti'^iUim does make the perturbation theory expansion of f/^^ 
difficult to see. Indeed, if we were to use Eq. ( |155| ) directly for this purpose, we would need to 
carefully keep track of the cancellations between contributions of the subtraction term and the 

contributions coming from the term \IAX~^t^DqIIA\ . For this reason we would like to find 



an alternative expression for t/^^ where all terms contribute with a positive sign. 

In order to expose the term in (u2~^r'^ Dq^^U^ that will cancel the subtraction term, we 
use the AGS equations, Eqs. (|126|) , to write 



(ui-'rfD^^u) = [(1 + UT)dfDoiI{TU + 1]] = 4m<i^oi 

\ / nm L J nm 

3 3 3 

+ ^ UnitidY^SimdfDoi + ^ d,fDoiSnktkdk^Ukm + ^ Unitid^%kdiDoitkd^^Ukm- (156) 

1=1 k=l l,k=l 

The last term in this equation can be written as 

3 3 

Unitid^^SikdfDoitkd^^Ukm = ^ Unitidi ^{5ii + 5ii)SikdfDoi{Sik + Sik)tkd^^Ukm- 

l,k=l l,k=l 

Expanding out the brackets we obtain four terms, three of which can be simphfied: 

3 

Unitid^^5iiSikdfDoi5iktkd^^Ukm = 0, (157) 

l,k=l 
3 

Unitid-i ^Siidikd^DQiSiktkdf^^iJkm = UniUd^ ^did^ ^ (jJim — GoSim^ , (158) 

l,k=l 

3 

^ Unltld^^iSlkdfDoiSiktkd^^Ukm = (Uni - GoSm^ d'^^d^d'^^tiUim- (159) 



Thus 



6 6 

(uX-htD-^U) - 5r,mdtD^i + V U^itid^^mdfDoi + V d'^DoiSnktkd^'Ukm 

+ UmUd-^dtd;^ (Uim - GoSim) + (t/ni - Go^i) di dj^ d^ ^iUim 
3 

+ ^nitidl^5ii6ikdlDQi5iktkd'f,^Ukm- (160) 

Using that 

3 

Unitidi ^5imd^DQi — ilnitid^ ^d^d^ ^GoSim 

1=1 

3 3 

= 5^ Ur,itidi%mdtD^,[l - 5u] = ^ Unitidr'6iJudfDoi (161) 

and similarly 

3 

di DoiSnktkd^^Ukm — GoSnid^ ^dfdi ^tiUim 

k=l 

3 3 

= diDoiSnktkd,^^Ukm[i — Sik] — d^ DoiSnk^iktkd^ ^Ukm, (162) 
fc=l fe=l 

we obtain that 




1=1 



k=l 



+2Unitidi ^d'^d^ ^Uim + ^ Unitidi ^Sii6ikd'^Doi6iktkdf^^Ukm- (163) 



Lk=l 



Substituting into Eq. ( |155| ) we finally obtain 



^ nm / J 



^nmdiDQi + Unltldi ^SimShd^DQi + d^ DQiSnk^iktkd^. ^Ukm 



k=l 



+ Uni {tiT^i + tfrfj ^) Ui^ + ^ ^SiiSikdfDoi5iktkd,^^Ukm 



Lk=l 



(164) 



This form for [/^^ has no subtraction terms and can be used directly to generate the cor- 
responding perturbation theory expansion. This is seen especially well from the graphical 
representation of Eq. ( |164|) given in Fig. |. 

Taking left and right residues of at the three-body bound-state poles leads to the three- 
body bound-state current (see Subsec. 5 below). In that case only the last three terms of 
Fig. H contribute. If we then consider the case of one-body currents with separable two-body 
interactions, we see that Fig. || reduces down to Fig. ||(a). 



4- j{ki) ijk transition current 

The strong interaction process j{ki) — > ijk where the final state consists of three free 
particles is described by the t-matrix 

Toj = Uoji/jj = UojDoj(j)j (165) 

where Uoj is given by 



To find the j{ki) ijk electromagnetic transition current Jqj (which can be used to describe 
processes hke pd — -ynpp), we proceed as before and define the Green function quantity 

It is Tqj which may now be gauged, thereby obtaining 

f,^ = G-,'%d-\ (168) 
Using the product rule we have that 




(169) 



and therefore the electromagnetic current can be written as 




5. Three-body bound-state current 

The three-body bound-state current was already discussed in Sec. Ill B above. There we 
obtained an expression, Eq. (0), which gives j'^ in terms of the two-body potentials Vi and 
the gauged two-body potentials . Here we would like to give an alternative expression that 
results from the gauging of the AGS equations. This has the advantage of giving j'^ in terms 
of the two-body t-matrices ti and the gauged two-body t-matrices . 

We recall that the AGS amplitudes Uij are defined through the expression for the 3^3 
Green function: 

G = GiSij + GiUijGj (171) 



where Gi is given by Eq. ( |114| ). Thus it is clear that Uij has a pole at the three-body bound 
state: 



where 



Xf = Gr^^p. (173) 



By writing Uij = Gq ^UijG^ ^ we may gauge Eq. (|1 71|) in the usual way. Then taking the left 
and right residues at the three-body bound state poles and using Eq. ( |137| ) we obtain that 

f = Y.^fG,T^,G,xl (1T4) 

Ik 

From Eq. ( |173|) it is easy to see that 



(175) 



where "^j" and are Faddeev components of the bound-state wave function as in Eq. (|87D , 
with ijk defined to be cychc. By introducing the cohimn matrix 



(176) 



with the same symbol \l/p being used as for the bound-state wave function, we may write 
Eq. ( p.74| ) in the matrix form 



;i77) 



giving us a formally identical expression to that of Eq. (|8T|) but where now each term on the 
RHS is a matrix. Interpreted as a matrix equation, this result expresses the current in terms 
of two-body t-matrices and gauged two-body t-matrices [see Eq. ( |140| )], while interpreted as a 
scalar equation (i.e. not a matrix equation) it expresses in terms of two-body potentials and 
gauged two-body potentials [see Eq. ([75|)]. 



6. {ijk) i{jk) transition current 

In the previous subsection we found the {ijk) {ijk) electromagnetic transition current 
by first gauging Eq. {T7l) for the green function G, and then taking left and right residues 



at the three-body bound-state poles. By contrast, in Subsec. 1 the j{ki) — > i{jk) transition 



current was found by first taking left and right residues of Eq. ( |1 71| ) at the two-body bound- 
state poles, which leads to Eq. ( |119|) , and then gauging this equation. It is straightforward to 
see that the final expressions for the electromagnetic transition currents do not depend on the 
order in which the gauging and the taking of residues is done. 

To determine the {ijk) i{jk) electromagnetic transition current it is convenient to first 
take the left residue of Eq. ( |171| ) at the two-body bound-state pole, then gauge the resulting 
expression, and finally take the right-hand residue at the three-body bound-state pole. Taking 
the left residue of Eq. ( |171| ) at the bound-state pole of particles j and k, but keeping the left 



propagator for particle i, leads to the Green function quantity 

Xij = (piGoUijGj = (piUijGQ'^Gj. (178) 

Gauging this equation, taking the residue at the three-body bound-state pole on the right, 
and then eliminating the left propagator of particle i, gives the {ijk) i{jk) electromagnetic 
transition current: 

= dr^ (^^fGoxf + '^^G.UuG.T^G.xdj ■ (179) 

This expression can be written using matrix notation as 

= ^f-qjp + $^ wr^^p (180) 

where $j and $f are row matrices whose j'th elements are defined by d^^(j)i5ij and d^^cp^Sij 
respectively. 



7. (ijk) — > ijk transition current 



The {ijk) ijk electromagnetic transition current jg describes the photodisintegration of 
the three-body bound state leading to three free particles, 'y{ijk) ijk. To find j^ we may 
start with the seven-point function of Eq. (0), take the right-hand residue at the three-body 
bound state, and multiply on the left by Gq^, in this way obtaining 

jI; = G.'GV^^p. (181) 

This gives the transition current for the photodisintegration process in terms of two-body 
potentials Vi and gauged two-body potentials ff . 

Alternatively, to obtain the corresponding expression in terms of two-body t-matrices ti and 
gauged two-body t-matrices , we may proceed as in the previous two subsections and start 
with Eq. ( |171| ) which may be written as 

G = SijGj + GiGq ^UijGQ ^Gj = SijGj + (1 + GoTi)UijGQ ^Gj. (182) 

Gauging the latter form of the equation, taking the right-hand residue at the three-body bound 
state and then multiplying on the left by Gg ^ gives 

= Gg 1 (G^T, + GoTt) Goxf + $^(1 + T,Go)UuTf,GoX^ ■ (183) 

Ik 

Now since (1 -|- TiG^jUn = Uqi, the electromagnetic transition current can also be written as 

jI; = (r^GoT, + Tt) Goxf + Yl Uoi^kGoX^ (184) 

Ik 

with Goxf given by Eq. ( |175| ). 



IV. SUMMARY 

In this article we have presented a general method for incorporating an external photon 
into a system of particles whose strong interactions are described nonperturbatively by integral 
equations. This method consists of gauging the integral equations themselves, and has the 
important feature of coupling the external photon to all possible places in the strong interaction 
model. As the photon is coupled everywhere, gauge invariance of all expressions for on-mass- 
shell electromagnetic transition currents is guaranteed. 

To discuss the details of our approach we have chosen the case of three distinguishable 
particles (with no coupling to two-particle channels) whose strong interactions are described 
by standard four- dimensional integral equations of quantum field theory. This type of three- 
particle system presents the simplest case for which no practical gauging method has so far 
been available. In the two-particle sector there have been previous gauging procedures that 
established the conserved currents for the NN system [|1| (no coupling to one-particle channels) 
and the nN system (coupling to the one-nucleon channel included), yet even in these cases 
the gauging of equations method provides a much simpler way to derive the same results. 

By gauging the integral equation for the three-body Green function where the kernel consists 
of two-body potentials Vi, we obtained an expression, Eq. (|115|), that describes all possible 



electromagnetic transition amplitudes of the three-body system in terms of the Vi and the 
gauged potentials v^. We have also shown how our method can be used to gauge the Alt- 
Grassberger-Sandhas equations for three particles in order to get more practical relations where 



the electromagnetic transition amplitudes are expressed in terms of two-body t-matrices ti and 
gauged t-matrices , see Sec. Ill C. 

Although we have presented the gauging of equations method within the context of quantum 
field theory where the integral equations are four-dimensional, it should be noted that the 
method itself can be used in a wider context. Indeed we have already used this method to 
incorporate an external photon into the three-dimensional equations of the spectator approach 



11,12|. Similarly, one could apply the gauging procedure to the three-dimensional approach 



of time ordered perturbation theory when the equations are expressed in terms of convolution 



integrals |22|. The gauging method also does not depend on the nature of the external field 
involved, so our results remain valid if the external field is due to a strongly or weakly interacting 
probe. In this sense the gauging of equations method provides the solution to the long standing 
problem of how to incorporate an external field into a nonperturbative description of quarks or 
hadrons. 
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